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P5 ' Abstract 

^3 : 

A perturbative study of a general class of lattice Dirac operators is reported, 
\ which is based on an algebraic realization of the Ginsparg- Wilson relation in the 

form 75(75-0) + (75.0)75 = 2cP'^~^^ ['^^D)'^^'^'^ where k stands for a non-negative in- 
^ . teger. The choice = corresponds to the commonly discussed Ginsparg- Wilson 

!>• \ relation and thus to the overlap operator. We study one- loop fermion contributions 

to the self-energy of the gauge field, which are related to the fermion contributions 
to the one-loop (3 function and to the Weyl anomaly. We first explicitly demonstrate 
that the Ward identity is satisfied by the self-energy tensor. By performing careful 
analyses, we then obtain the correct self-energy tensor free of infra-red divergences, 
• as a general consideration of the Weyl anomaly indicates. This demonstrates that 

our general operators give correct chiral and Weyl anomalies. In general, however, 
the Wilsonian effective action, which is supposed to be free of infra-red compli- 
p • cations, is expected to be essential in the analyses of our general class of Dirac 

D . operators for dynamical gauge field. 
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1 Introduction 

Recent developments in the treatment of fermions in lattice gauge theory are based on a 
hermitian lattice Dirac operator 75D which satisfies the Ginsparg- Wilson relation[l] 

75/^ + 1^75 = 2aD^^D (1.1) 

where the lattice spacing a is utilized to make a dimensional consideration transparent, 
and 75 is a hermitian chiral Dirac matrix. An explicit example of the operator satisfying 
(1.1) and free of species doubling has been given by Neuberger[2]. The relation (1.1) led 
to an interesting analysis of the notion of index in lattice gauge theory[3]. This index 
theorem in turn led to a new form of chiral symmetry, and the chiral anomaly is obtained 
as a non-trivial Jacobian factor under this modified chiral transformation [4]. This chiral 
Jacobian is regarded as a lattice realization of that in the continuum path integral [5]. See 
Refs.[6] for reviews of these developments. 
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We have recently studied a specific generalization of the algebra (1-1) [7] 

75(75^) + (75^^)75 = 2a''+\j,D)''+' (1.2) 

where k stands for a non-negative integer and k — corresponds to the ordinary Ginsparg- 
Wilson relation. When one defines 

H = ^5aD (1.3) 

(1.2) is rewritten as 

^,H + H^, = 2H^^+\ (1.4) 
The algebra (1.4) is equivalent to a set of equations 

//'75 - 75/^' = 0. (1.5) 

where the second relation is shown by using the defining relation (1.4), and the first 
relation in (1.5) becomes identical to the ordinary Ginsparg- Wilson relation (1.1) if one 
defines H(2k+i) — H'^^'^^- One can thus construct a solution to (1.5) by following the 
prescription used by Neuberger[2] 

where 

Df"-"-^ = + 3'"'+' - (I^)2fc+i. (1.7) 

Here wc note that the conventional Wilson fermion operator Dw ( with a non-zero mass 
term ) is given by 

Dw{x, y) = i^Ci^ix, y) + B{x, y) - -mo5:c,y, 

a 

U^{y) = exp[iagA^{y)]. (1.8) 

The parameter r stands for the Wilson parameter. Our matrix convention is that 7^^ are 
anti-hermitian, (7'*)^ = —7^, and thus = 7'*C^(n, m) is hermitian 



(1.9) 



The hermitian operator H itself is then finally defined by (in the representation where 
H{2k+i) is diagonal) 



H ^ {H^,,+uy/''+' (1.10) 
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in such a manner that the second relation of (1.5) is satisfied, which is in fact confirmed 
in the representation where H(2k+i) is diagonal[7]. Also the conditions < mo < 2r — 2 
and 

2mf+'' = 1 (1.11) 

ensure the absence of species doublers and a proper normalization of the Dirac operator 
H. 

The locality properties are crucial in any construction of lattice Dirac operator, and the 
locality of the standard overlap operator with A; = has been established by Hernandez, 
Jansen and Liischer[10], and by Neuberger[ll]. 

As for the direct proof of locality of the operator D for general k, it is shown for 
the vanishing gauge field by using the explicit solution for the operator H in momentum 
representation[12] [9] 

H{ap^) = 75(^)^(^)^{(^/^ + M,)^-(y^-M,)^^} 

, 1 . _fc+l_ . 1 . _fc±L r / / ~ . fc + 1 



J. fc+1 X fc+1 / fc+1 / " fc 

= 75(:^)^(^=)^{(V^(fe) + Mfe)^ - ( -Mk)^ 



;i.l2) 



where 



and 



S/j, = sm apij, 
Cjj, = cos apij, 

^=7^sinap^. (1-14) 

For k = 0, this operator is reduced to Neuberger's overlap operator [2]. Here the inner 
product is defined to be > 0. This operator for an infinitesimal p^, i.e., for \ap^\ <C 1, 
gives rise to 

H ~ -75a x<l + 0(ap)') +75(750 (1.15) 

to be consistent with H — 750!); the last term in the right-hand side is the leading term 
of chiral symmetry breaking terms. The locality of this explicit construction (1.12) has 
been shown by examining the analytic properties in the Brillouin zone[12]. 

It is important to recognize that this operator is not ultra-local but exponentially 
local[13]; the operator H{x,y) in (1.12) decays exponentially for large separation in coor- 
dinate representation as [12] 

H{x,y) ~ exp[-|x-7/|/(2.5/ca)]. (1.16) 
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An explicit analysis of the locality of the operator H{2k+i) = H'^^^^ {not H itself ) in 
the presence of gauge field, in particular, the locality domain for the gauge field strength 
has been performed. The locality domain for ||-F)tjy|| becomes smaller for larger /c, 
but a definite non-zero domain has been established[12]. The remaing task is to show the 
locality of the operator H = {H(2k+i)Y^^'^^^^'' itself in the presence of gauge field. Due to 
the operation of taking the {2k + l)th root, an explicit analysis has not been performed 
yet, though a supporting argument has been given in Ref. [12]. 

It has been shown that all the good chiral properties of the overlap operator[2] are 
retained in the generalization in (1.4) [8] [9]. The practical applications of this generalized 
operator D are not known at this moment. We however mention the characteristic prop- 
erties of this generalization: The spectrum near the continuum configuration is closer to 
that of continuum theory and the chiral symmetry breaking terms become more irrelevant 
in the continuum limit for /c > 0. The operator however spreads over more lattice points 
for larger A;, as is indicated in (1.16). 

In this paper we study a perturbative aspect of the general class of Dirac operators. 
To be specific, we study the one-loop fermion contribution to the gauge field self-energy, 
which is related to the (3 function and to the Weyl anomaly. 



2 Self-energy tensor, /^-function and Weyl anomaly 

The lattice perturbation theory is very tedious in general[14]-[22], and it is more so in our 
generahzation. For this reason, we study the simplest diagrams related to the one-loop 
self-energy correction to gauge fields. This effect is also related to the fermion contribution 
to the lowest order /3-function and to the Weyl anomaly [23] [24]. A rather general analysis 
of Weyl anomaly is possible, and we first briefly summarize it. 

In the standard continuum formulation, one starts with the path integral defined in a 
background curved space [25] 

j dn&cp[j d^x^i^iiPi^]. (2.1) 
The general coordinate invariant path integral measure is defined by 

d^l = V^Vi) (2.2) 

and the Weyl transformation laws are given by 

]^a{x)\i^{x), 

\a{xmx) (2.3) 

where e'^{x) stands for the vierbein. This transformation law is fixed by the invariance of 
the action in the above path integral under a global (i.e., constant) a, and the Weyl weight 
factor of fermionic variables is essentially defined by the vierbein in ^ = e(^(x)7"D^. 



e^(x) ^exp[Q;(x)]e^(x), 
i^{x) = {gy^'^^|J{x) exp[- 

i'{x) = {9f''^'4^{x) exp[- 
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The Jacobian for this transformation of fermionic variables is given by 
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In J(q;)= lim Tra(x)exp[-(p/Mf]^Tra(x)-^F^^F'"' (2.4) 

M— >C» 24:71 

where the mode cut-off of p is provided by e"^'^"/^-''^ in terms of the eigenvalues of p. 
See Ref.[25] for further details. 

When one analyzes the higher derivative theory 

Ait+i = / d^x^i^iiPf'+'Tj;, (2.5) 

the Weyl transformation laws are given by 

e^{x) exp[a{x)]e^{x), 

2/c + 1 ~ 
i'ix) = {gY^'^ipix) exp[ — a{x)]ilj{x), 

~ _ 2A- -I- 1 ^ 

i^{x) = {gy/'^|^{x) ^ exp[ ^aix^x) (2.6) 

and the Weyl anomaly is given by 

In J2k+i(a) = lim Tr(2k + l)a(x) eM-iiP?"^' /M^'+^n (2.7) 

M— >oo 

Since the Weyl anomaly is independent of the regulator function [25] , we have 

lnJ2fc+i(Q;) = (2A; + 1) In J(q;). (2.8) 

This relation (2.8) is also understood from a view point of the self-energy correction 
to the gauge field as follows: 

det(^)2'=+^ = exp[(2A; + l)Trln^] 

= exp[(2A; + l)Trln(^-i^/)] (2.9) 

= exp[(2A; + l)rrln ^ ^{igfTr- 4^4 + ....] 



The term quadratic in the gauge field Aj^ gives the self-energy correction, which is 2A; -|- 1 
times larger than the self-energy correction generated by det p. 

This analysis of the self-energy correction is applicable to the present lattice operator. 
By our definition in (1.10) we have 

exp[Trln//] = exp [Tr In //(VSr^l = ^M^j^Tr In H^^k+i)]- (2.10) 
For a sufficiently small coupling constant g, we have 

exp[Tr In H] = exp[Tr In + Trg^A^{x)Oix, y)^.A,iy) + Oig^)] (2.11) 
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where H^^^ stands for the free Dirac operator given in (1.12), and the second term stands 
for the lowest order term in the effective potential and thus for the lowest order self-energy 
correction to the gauge field. Similarly, we have 

exp[^^Trlni/(2fe+i)] (2.12) 

^ ^^p[^^^^^"^Sh-i) + ^^Trg'A,(x)d(x,y),^A^(y) + 0(g% 

where Tr In H^2k+i) stands for the free part of H(^2k+i) ■ Those zeroth order terms satisfy 
the relation 

Tr In (°) = 2^Tr In H^^l,^ (2.13) 

if one uses the explicit form of the operator in (1.12). 
We thus conclude 

Trg^A^{x)0{x, y)^.A,{y) = ^^Trg^A^{x)d{x, y)^.A,{y) (2.14) 

for a sufficiently small coupling constant g, which shows that the lowest order self-energy 
correction in the left-hand side for the operator H is evaluated by the self-energy correction 
in terms of H(2k+i)- Wc use this relation for the evaluation of the lowest order self- 
energy correction for any k >1. Note that the operator i?(2fc+i) = H'^^^^ is much better 
understood than H itself in our construction. We also confirm that this relation (2.14) is 
in fact valid by evaluating the left-hand side directly for the simplest case k — 1. 
^From a view point of Weyl anomaly, one may tentatively take 

lim Tri2k + l)a{x) exp[-(i7(2fe+i)/(aM)2'=+^)2] 

lim Tr(2A;+ l)a(a;)exp[-((^)2'=+VM2^+^)2] (2.15) 

M— >oo 

as a lattice version of the Weyl anomaly. We then obtain the same result as the self-energy 
correction in the limit a — > 0, although no systematic formulation of Weyl anomaly on 
the lattice is known. 

In lattice perturbative calculations, however, we should be careful of the possible 
appearance of infrared divergences, which should cancel in the final result. We show that 
a careful analysis gives the correct result of continuum theory free of infrared divergences 
for a — > 0. 



3 The vacuum polarization tensor by H(^2k+i) 

In this section, we calculate the one-loop fcrmion contribution to the vacuum polarization 
U^i, on the basis of the operator H(2k+i) ( and not H itself) following the analyses in 
Section 2. We first show that the Ward identity is satisfied to be consistent with gauge 
invariance and that there appear no divergences except for the logarithmic divergence for 
a — > 0. We then discuss the gauge field wave function renormalization factor. 
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Feynman diagrams for the vacuum polarization witli fermion one-loop are shown in 
Fig.l, and the necessary Feynman rules are given in Appendix A. 

The amplitude corresponding to Fig. 1(a) is given in terms of the notation in Appendix 
A by (by using tr{T'^T^) = \5^^ and Nf flavors in QCD) 



nS)(p) 



-9^ Nf.AB 



xtr 



4a4'=+2 2 Jt{w{t)+w{t + p)y 



D^\th, |xi^(t, t + p, -p) - ^Xlit, t + p, -p) 



xD^,\t+p)jAXi,{t + p,t,p) 



w{t) 
Xo{t + p) 
w{t + p) 



Xl{t + p,t,p) 



Xo{t + p) 
w{t + p) 



w{t) 



We omit the factor S^^ from now on. 

The amplitude corresponding to Fig. 1(b) is similarly given by 
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202^=+! 2 Jt 



Ni f 
2 Jt 



tr 



2w{t) 



X 



X2,At,t, -p,p) - ^xUt,t, -p,pf'^^^ 



w{t) 



w{t) 



(3.1) 



{w{t) + w{t + p)y 

X {Xi^{t,p + t, -p)Xl{p + t, t,p)Xo{t) 
+Xi^(t,p + t, -p)Xl{p + t)Xi,(p + t, t,p) 
+X^{t)Xl{t,p + t, -p)XUp + t, t,p) 
2w{t) + w{p + t) 

w{tyw{p + 1) 

xXo{t)Xl{t,p + t, -p)Xo{p + t)Xl{p + t, t,p)X^{t)} 
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X it] X (t) 

w[t) ^ w[t) 

1 



{wit) + w{t-p)y 

X {Xi,(t, t - p, p)Xl{t - p, t, -p)Xo{t) 

+Xi,(t, t - p, p)Xl{t - p)X^^{t - p, t, -p) 
+Xo(t)Xj,(t, t - p, p)Xi^{t - p, t, -p) 
2w{t) + w{t-p) 



w{tyw(t — p) 



xXoit)Xl{t, t - p,p)Xoit - p)xUt - p, t, -p)Xo 



(3.2) 



3.1 Ward identity 



We first show that the Ward identity for 11^,^ as a manifestation of gauge invariance holds 
as follows, 

Ep^(n!rj(p) + n!f)(p))) = 0, (3.3) 

V 

where py = (2/a) sinap;^/2.[| For this purpose we first calculate J^uPuXii, and J2uPuX2^u- 
For J^^PuXii, we have 



l+m=2k 

+ (^Y.il-cos{t + p)pa)^ |(^r5](l-cos(t + p)pa)^ - ^(1 - cost^a)^ | 
X r'^{\ — costptt) 

I \ 2fc+l 

= lii^^pf""^^ + [r E(l - cos(t + p),a) j - mf^^ 

{/ \ 2fc+l 

= a2'=+M^o(t + p)-Xo(t)}, (3.4) 
and further we have 

Y.p^ fx,.(t + p,t,p) - ^^^xL(t + p,t,p)M)) 

V V wit^p) wit) j 

= 275a^'+^(w(t + p) + u;(t))(/^o(p + t) - ^o(t)), (3.5) 

^The Ward identity in the case of the overlap Dirac operator has been confirmed exphcitly in Ref . pO| . 



where we used the following relations 



^p^7^ cos(i + p/2)^a = -{^t+p - h 

1 1 Oj 



^Pj,rsin(t + p/2)^a = ^ I ^r^(l - cos(t + p)pa)j - ^(1 - costpo) 



Using the above relations, we also have 



E 

l+m=2k L 



iiih)^ ( cos ( t + I ) a ) (i^+p)' 



r sin ( t + ^ ) a ) ( r 



^(1 -cos(i+p)pa) 
p J 



-i{i^t)^ |^i7^cos (^t + aj {i^t 
+ ^r^(l-costpa)j 1^ 
— ^'^XK-'- ~ cosipo)^ ^rsin + ^ aj ^(1 — costpa) 
+ ^ 1^ (^/t)' (^^7m cos + I) « - ^7// cos (^t - (V*)"" 
+ ^r^(l - costpa) j + a 

X ^r^(l -cosipa)j | 



m 



— r sin ( t — - ) a 



Therefore we obtain 



E [^2/.!/ (t, t , -P, P) + X2ut, {t, t, p, -p)] 

V 

^ Xi^{t,t + p,-p) - Xii^{t-p,t,-p). 



By using these relations, Pi^^l^^ (p) is written as 

-g^ Nf r 1 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



E^^n(t)(p) 



2a2*^+i 2 Jt {w{t) + w{t + p)} 



■tr 



{D^\t)-D^\t + p)}^, 



><{^,^(,,^,,_,)_^^t^(,,^,,_,)^ 



(3.10) 
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and similarly Y^uPi^^^^Kp) written as 



tr 



X 



1 



w{t) + w{t — p) 



- p, t, -p) - -^;;^Xl{t - p, t, -p) ^ I 



1 



1 



^tr [{D„-i(t)-Do-i(t + p)} 



75 



X 



w(t) + w(t + p) 



Xi^{t, t + p, -p) - ^^xlit, t + p, -p) 



Xo{t + p) ' 
w(t + p) 



(3.11) 



Combining these two expressions, the Ward identity for the vacuum polarization tensor 
holds as in Eq. (|3.3|) . This Ward identity dictates the tensor structure of n^,^(p) for small 
P^ to be 



^,j,u{p) ^ {p 6f,^ - Pi,pt,)Il{a p ). 



(3.12) 



3.2 Structure of divergences 

We next examine the structure of various divergences. To evaluate the divergent parts 
of n^")(p) and n|fj(p), we rescale the integration momenta t^/a in each amplitude 
( ^ and (|3]|). For QCD with Nf flavors, we obtain 



n!rj(p) 



1 



tr 



w{t) + M{t) 
X <^ Xi^{t,t + pa, -pa) 



t {wit) + w{t + pa)} 
Xoit) 



\w{t + pa) + M{t + pa) 
X < Xi^(t + pa,t,pa) 



wit) ^^'^ 
1 



XL{t,t + pa,-pa) 



Xo(t + pa) 
w{t + pa) 



—— -Xi^{t + pa,t,pa)—— 

w[t + pa) w[t) 



(3.13) 



and 



n?J(p) 



Aa"" 



tr 



w{t) + M{t) 



+ 1 



2w{t) 
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X 



X2^,u{t,t, -pa, pa) - 
1 



w{t) 



xLu{t,t, -pa, pa) 



Xojt) 
wit) 



{w{t) + w{t + pa))'^ 
X |Xi^(t,pa + t, -pa)Xl^{pa + t,t,pa)Xo{t) 

+Xi^{t,pa + t, -pa)Xl{pa + t)Xiy{pa + t,t,pa) 
+Xo{t)Xl^{t,pa + t, -pa)Xi^{pa + t,t,pa) 
2w{t) + w{pa + t) 
w(t)'^w{pa + 1) 

xXo{t)Xl^{t,pa + t, -pa)Xo{pa + t)Xl^{pa + t, t,pa)Xo(t)| 

X it) X it) 

+X2uf, (t, t,pa, -pa) - -^^Xl^^ (t, t,pa, -pa) - ° 



w{t) 



wit) 



{w{t) + w{t — pa))'^ 
X |Xij,(t, t — pa,pa)Xl^{t — pa, t, — pa)Xo(t) 

+Xii,{t,t - pa,pa)Xl{t - pa)Xifj_(t - pa,t, -pa) 
+Xo(t)Xl^{t, t - pa,pa)Xif,{t - pa, t, -pa) 
2w{t) +w{t-pa) 
w(t)'^w(t — pa) 

xXo{t)Xl^{t,t - pa,pa)Xo{t - pa)Xl^{t - pa,t, -pa)Xo{t)} 



(3.14) 



where = J'^^dH/{2n)'^. In the above two equations w, Xq, Xi^, X2^ij are appropriately 
redefined according to the rescahng of t^. For example, 



w{t) 



\ 



(sf) 



2\2fc+l 



+ 1 (^E(l-cost, 



2fc+l 



2fc+l 



(3.15) 



where sf = I]/^ sin^ 

We first want to show that there are no nonlocal divergences of the forms p^ /{a^p^)"^ 
or p^Pv/ {a^p^)^ {n > 2). For this purpose we confirm that Eq. (|3.13| ) and Eq.( |3.14D are 
not singular for p = 0. Setting p = in these equations, we have 



n!rj(o) 



1 



tr 



w{t) + M{t) 



4w(t)2 



x{x,(M.O)-|fx.(M,0)^ 
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and 



\w{t) + M{t) 
X |xi,(t,t,0) - 



—^Y^Xl^{t, t, 0)—4-r- 
w[t) w{t) 



nS(o) 



4a2 



Hi 



w{t) + M{t) 



+ 1 



2w{t) 



X 



t, 0, 0) - ^xl^t, t, 0, 0)^ 

w[t) w[t) 



X {Xi^{t,t,0)xl{t,t,0)Xoit) 

+Xi^{t,t,0)X^oit)Xiu{t,t,0) 

+Xo{t)xl{t,t,0)X^,{t,t,0) 
3_ 

~w{ty 

xXo{t)Xl{t, t, 0)Xo{t)Xl{t, t, 0)Xo{t)} 



1 



Xojt) 
w{t) 



xl^{t,t,o,o) 



Xojt) 
w{t) 



x{x^,{t,t,Q)Xl{t,t,0)Xo{t) 

+X^,{t,t,0)Xl{t)X^^{t,t,0) 

+Xo{t)Xl{t,t,0)X^^{t,t,0) 
3 

xXo{t)Xl{t, t, 0)Xo{t)Xl{t, t, 0)Xo{t)} 



(3.16) 



(3.17) 



Now on the basis of the expressions of w, Xq, Xi^ and X2fj,u and the fact that there are 
no doublers, the possible singularity may occur only around the region t ~ in each 
integral. Only the fermion propagators can exhibit singular behavior for t ~ 0. The 
leading singularity in t in n^"^(0) vanishes as 



n^tHo) 



tr 



it') 



k+l 



it') 



k+l 



(5<1), 



(3.18) 
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and similarly the leading singularity in n|f2 (0) vanishes as 



nS(o) 



tr 



2\k+l 







(3.19) 



Higher order terms in t are obviously non-singular. Since both n|j")(0) and 11^^2(0) are 
not singular, Ii^y{p) does not have the non-local divergences of the forms |}^/(a^p^)" or 

There may still exist the quadratic divergence in I\.^y{p). From Eq.( p.3| ), the form of 
the quadratic divergence for small is 



(3.20) 



with a constant C. We have already established that the ^ limit of Il^^{p) is well- 
defined, which excludes the singular term p^Pu/p^] this term depends on the direction of 
the approach — > 0. We thus conclude C = 0, namely, the quadratic divergences cancel 
between diagrams (a) and (b). 

Next we confirm that there are no divergences of the structure such as a^p"^ x oo, etc, 
which vanish in the naive continuum limit. These unusual divergences , which may be 
termed as infrared singularities, may occur in our treatment of H{2k+i) which corresponds 
to a higher derivative theory on the lattice. These divergences, if they should exist, could 
appear in the integration region around t ~ and could remain even for arbitrarily small 
p. Therefore we evaluate Iljfj (p) (g]T^) and ng](p)(^]T|) for < 6^ and ap ~ 0. After a 
straightforward calculation, we obtain 



n!rj(p) 



2a2 



tr 



1 



\2k+l 



A+m=2k 



nS(p) 



Nfg' 



a+m=2k 



(3.21) 



2a2 7*2 <52 



tr 



\2k+l 



l+m+n=2k~l 



l+m+n=2k-l 



(3.22) 



These amplitudes (a) and (b) separately could contain infrared singularities. The cancel- 
lation between the amplitudes (a) and (b) further takes place as: 

n?j(p) + nW(p) ^ 



2a2 



-(2A; + 1) / tr 



1 . 1 



(3.23) 



This final expression, which has the same structure as that in continuum theory, means 
that there are no divergences such as a^p^ x oo, etc. 
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Finally, we investigate the logarithmic divergence. From the above analyses, we know 
that Ilfj^iy{p) does not have the divergences of the negative power in a. Therefore if there is 
the logarithmic divergence in 11^,^ (p), it appears from the singular part in the integral for 
a — ^ and thus the singular part should appear in the integration region around t ~ 0. 
We first evaluate U^^J (p) . There are several ways to extract the logarithmic divergence [B^ 



22 1 . Here we use the procedure discussed in the paper by Karsten and Smit [21]. First 



the denominators of the propagator are combined using Feynman parameters and the 
integration variables are shifted t^— p^ax as follows. 



~ 8a2 (r(2A; + l))' l-.+pax (2^ io ^ {w{t - pax) + w{t + pa{l - x))}- 



X 



[a{l - x) + (3xf^^^ 



xtr 



si 



■pax) 



't—pax 



pax , t + pa{l — x) , —pa) 

Xo{t-pax) J _ , , _ ^^ . Xo(t + pa(l 
w{t - pax) ^'^^^ pax,t + Ml P^K{t + pa{l 

^{slpa{l-x)f h+pa(l-x) {Xi^{t + pa{l-x),t 

Xo{t + pa{l - x)) 




w(t + pa{l — X 



pax, pa) 

X)) , ,^ , .Xoit — pax)] 

—^Xl{t + pail — x),t — pax, pa) — — > 

;)) w{t — pax) J 



where^ 



a = {w{t — pax) + M(t — pax)}^^^'^''^^\ 

(3 = {w{t + pa{l-x)) + M{t + pa{l-x))}^^''^''^^^ 

Then we split the integration domain into two regions as follows 



+ 



5< 1, 



(3.24) 



(3.25) 
(3.26) 



(3.27) 



and we evaluate the t"^ < 6^ part in the continuum limit, ignoring the t"^ > 5"^ part which 



does not contain divergence. Eq.( |3.24D is the complicated integral including sines and 
cosines. However for < 5^ and a ^ with fixed small we can expand both the 
denominator and the numerator of Eg. ( p. 24 ) separately in powers of t and a, and we have 



nS(p) 



dx 



X 



{1-x 



-Nfg"^ r(4A; + 2) r dH 
2a2 (r(2A; + l))' lt^<5'^ (2^^ 7^ + ^2^23,(1 _ 3,)} 



I 4fc+2 



^We used the Feynman's formula: 

1 _ r(4fc + 2) 



a 



2fc+l ^2fe+l 



r(2fc + i)r(2/c + i) 



dx 

\ax 



x'^{l-x) 



2k 



+ /3(l-x)] 



4fe+2 • 
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xtr 



0<l<2k 



[o<m<2fc 

(3.28) 

We next evaluate n|f2(p) in the similar way and we obtain 

-Njg^ T{4k + 2) f dH x'^^{l-xf^ 



^^^^ - 2a2 (r(2fc + 1))' U<5^ (27r)4 k {t^ + p^aH{l - x)f''+^ 
xtr + i)a{l - x))f^^^ I - i)ax)f^^^ 



X { E iW -i)ax)\ i-t^ 

[o<Z+m<2fc-l 

X + i>a{l - x))r ^7. W - i>ax)f''-'-'-"'} 



X 



E i +|^a(l - a;))]'i7t. 

0<Z+m<2fc-l 



By this way n^j,(p) = Ii^^}{p) + n^''2(p) written as 

-Nfg^ TUk + 2) r dH , x^'^fl - xf' 

' ' dx 



(3.29) 



2a2 (r(2A; + 1))' U<s^ {2n)^ Jo {^2 + ^2^2^(1 _ ^)}4fe+2 
X {2k +l){t- paxY''{t + pail - x)YHr [(/ - ^ax)-f^{f + ^(1 - x))-f^] (3.30) 

The singular part corresponding to the logarithmic divergence is obtained from the leading 
part in t and a. Noting the spherical symmetry of the integral and dropping 0{a^) terms 
in the numerator, the singular part is given by 

^ , , -Nfg^ f dH , x'^Nl-xf' 



2a2 A2<52 (27r)4 Jo {^2 + ^2^2^(1 _ 
x4(2A; + l) [-1(^2)4'+ V 

'4.. 4 



+ (t2) Va^V {-(^' + (1 - ^)') {^k' + ^A;) + a:(l - x) (^fc^ + 2fc + l) } 
+(i')'W^«' {(^' + (1 - ^)') (^^' + - ^(1 - ^) {Ik' + 4A; + 2) 



(3.31) 
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After some calculations, the term proportional to logp^a^ is obtained as (by restoring the 
factor S^^) 

{2k + l)6^^^{p\^ - p^p^) \ogp'a\ (3.32) 

Combined with the general analysis (2.14) in the previous section, we conclude that the 
divergent part of the gauge field renormalization factor arising from fermion one-loop 
diagrams for the general Dirac operator D — {■^^/ a)H is given by 

^A=l + ^log/xV, (3.33) 

where is the renormalization scale. This factor indeed reproduces the correct result for 
the QCD-type continuum theory [23] [24]. 

Incidentally, the result (3.32) could also be directly obtained from (3.23), which cor- 
responds to 2A; -h 1 times the vacuum polarization tensor generated by a conventional 
massless fermion. 



4 The vacuum polarization tensor for H with k = 1 



In this section we calculate the one-loop fermion contribution to the vacuum polarization 
tensor 11^,^ on the basis of H with the simplest case k — 1. We perform essentially the 
same analysis as in the previous section. 

Feynman diagrams for the vacuum polarization with a fermion loop are shown in Fig.l, 
and the necessary Feynman rules are given in Appendix B. 

The amplitude corresponding to Fig. 1(a) is given by (for QCD with Nf flavors) 



■AB 



/■ 

Jt I 



t a{t, t + pY 
D^\t) {Z)(t, t + p)75//(3)i^(t, t + p, -p) 

-75Ho{t){75Hi3)i^,it,t + p, -p))h5Ho{t + p)} 
xD^\t + p) 

X {D{t + p,t)-f5H(^3)i^{t + p,t,p) 

-75//o(^ +P)(75^(3)i.(^ +P,^,p))^75^o(i)} 



The amplitude corresponding to Fig. 1(b) is 



tr 



1 



a{t, t) 



X {D{t,t)'yr^H(^3)2f,u{t,t, -p,p) 
-75-f^o(^)(75-H'(3)2Mi^(^, t, -p,p))hbHQ{t) 



(4.1) 
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-Hoity^^H.^it, t + p, -p){^,H,,{t + p, t, p))h5Ho{t) 

-2Ho{tf^,H,^{t, t + p, -p){^,Ho{t + p))^,Hi,{t + p, t, p) 

-Ho{tf-f5Ho{t){-f,H^^{t, t + p, -p))^^^Hi,{t + p, t, p) 

+^5Ho{t)i-f5Hi^{t, t + p, -p))t75i/o(t + p)il5HUt + p, t, p)y-f5Ho{t) 

+ ip,fx^ -p,u)}], (4.2) 



where 



-f5Hi^{t+p,t,p) = -^^^-p-y {L'(t + p,t)75if(3)i^(t+p,t,p) 

-IbHoit + p) (75^(3)ii.(i + p, t, p))'^i^Hq 

(4.3) 

We first show that the Ward identity for n^^, holds in this case also. /^From the analysis 
in Subsection 3.1, we obtain 

^ Pu {l5H(3)2^^ {t, t, -p, p) + 75^^(3)21.^. (t , -P) 

V 

= 75^^(3)1/. (^, ^ + P, -P) - 75^^(3)1/. (^ - P, t, -p). 



Using these relations Y.uPu^^nl{p) is written as 



J, (aj ^ J t 



1 



a{t, t) 

X {D(t,t)(75i/(3)i^(t,t + p, -p) -75ff(3)i^(t-p,t, -p) 
-75ifo(^)(75^^(3)lAt(i, i + P, -P) - 75^^(3)1m(^ " P> ^> -p))"^75^^o(i) 

+Ho{tf-f5HiMt + P)Ho{t) - / X 75ifi^ 
-2Ho{tf^5Ho{t)HlHoit + p) + D(t, t + p)-iMt)HiMt) 
-H,{tf^,H[^H,{t-p)H,{t) +I'x ^,H[^ 

+2Ho{tfirMt)HlHo{t - p) - D(t, t - p)^Mt)H[^Ho 
where '/' means p — > — p and 



X {2/7o(t)'(^^(3)iM(^, ^ + -P) - ^(3)iM(i - p. t, -p) 
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(4.4) 



(4.5) 



I = -3Hoit)^ + 2HoitfD{t,t + p). (4.6) 
Noting Eg. ( |B.5| ), Y^uPu^I^IKp) is rewritten as follows, 



-Ho{t){H(^sMt, t + p, -p) - Hf^sMt - p, t, -p))Ho{t) 
+?,Ho{tfH^^{t, t + p, -p) - 2HoitfH^s-j,^{t, t + p, -p) 
+Ho{t)H^:i)i^{t,t + p, -p)Ho{t) 

-3Ho{t)^H,^{t - p, t, ~p) + 2H^{tfH(^,)^^{t - p, t, ~p) 
-Ho{t)H(^3)it.{t-p,t,-p)Ho{t)}' 

tr \D-\t)-f,{H,^{t, t + p, -p) - Hi^{t - p, t, -p)) 



a 2 



(4.7) 



We next calculate Y^u Pu^^/^J {p) , 



V 

X 



tr [Do\t)^5Hi^,{t, t + p, -p)D^\t + ph5{Ho{t + p) - Ho{t)) 



a 2 
a 2 



AB 



AB 



iD,\t) - D^\t + p))^,H^^{t,t+p, -p) 
D^\th,{H,^{t, t + p, ~p) - Hi^it - p, t, -p)) 



/^From Eg. ( [4.8|) and Eq.( [4.7|) , one can see that the Ward identity for the vacuum polar- 
ization tensor holds. 

We next examine the structure of various divergences. Rescaling the integration mo- 
menta t^/a in each amplitude, we obtain (by omitting the factor 6^^ from now 
on) 



nS(p) 



xtr 



2o? Jt a{t, t + pay 
Do ^(t) \^D{t, t + pa)75if(3)i^(t, t + pa, -pa) 

-75if (t) (75^{3)iM {t,t + pa, -pa) )^-f5Ho{t + pa)} 
xDQ^{t + pa) 

X |D(t + pa, t)75-f^(3)ii.(i + pa, t, pa) 

--f^Hoit + pa)(75if(3)i^(t + pa, t,pa))^-f5Ho{t)} 



(4.9) 



and 



nS(p) 



2a2 



tr 



Do\t) 



a{t, t) 



{D{t, t)'^^H(^)2i,v{t, t, -pa, pa) 
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-75-f^o(^)(75-H"(3)2Mi^(^, t, -pa,pa)y-f5Ho{t) 
-Ho{ty-f5Hiij_{t,t + pa, -pa){-f5Hi^{t + pa,t,pa))^-f5Ho{t) 
-2Ho{t)^-f5Hii^{t, t + pa, -pa){-f5Ho{t + pa)y-f5Hi^{t + pa, t,pa) 
-Ho{tf-f5Ho{t){'y5Hii^{t,t+pa, -pa))"^ -f^Hi^{t ^ pa,t,pa) 
+-f5Ho{t){-f5Hiu_{t, t + pa, -pa))^-f5Ho{t + pa){-f5Hi^{t + pa, t,pa)y-f5Ho{t) 
+{p,f,^-p,u)}]. (4.10) 

First, we want to show that n|^"j(p) and n|fj(p) are finite and well defined for p — 0, and 
thus the divergent terms of the forms and p^Pi,/ (a^p^)""^-*^ with n > do 

not appear. Setting p = in these expressions, we have 



nS(o) 



xtr 



-Nfl 
2a2 



a{t,ty 



-75i^o(t)(75%)i/.(t,t,0))S^o(t)} 

X {D{t,t)^5H^3-ji,{t,t,0) 

-75i/o(i)(75^^'(3)i.(i, 0))^75//o(i)} 



(4.11) 



and 



niS(o) 



2a2 



tr 



a{t, t) 

x{D(t,t)75i/(3)2^,(t,t,0,0) 
-75i/o(t) (75^^(3)2;.. (t, t, 0, 0))hMt) 

-Ho{tf^5Hi^{t, t, 0)(75i/i.(t, t, 0))h,Ho{t) 
-2H,{tf^^H,^{t, t, 0){J,Hoit)y^,H^,{t, t, 0) 
-Ho{tf-i^Ho{t){-i^Hi^{t, t, 0)y-f,H,,{t, t, 0) 
+j,Ho{t){j,H,^{t,t,0)y^,Ho{t){^5H,,{t,t,0)y^,Ho{t) 



(4.12) 



The singularity may occur around the region t ~ in each integral. However, the leading 
order part in t in 11^^") (0) and Il^^j (0) vanish as 



n?2(o) 



2a2 Jfi<&2 (t2)4 

Nfg' 



X tr 



0, 



~ 0, 
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for 5 <^ 1. Since both ll^^^{0) and H^^'jlO) non-singular, U^^ does not contain the 
non-local divergences. /^From the fact that the ^ limit of 11^,^ (p) is well-defined 
and finite and that the Ward identity holds, we also conclude that the possible quadratic 
divergences cancel between diagrams (a) and (b). See also the analysis in Section 3. 

Next we confirm that there are no divergences of the structure such as a^p^ x cxd, 
etc, which vanish in the naive continuum limit even if they existed. For this purpose we 
evaluate 11^") (p) and Il^^}{p) for < 5^ and ap ~ 0. We thus examine the behaviour of 
various functions appearing in these amplitudes for t"^ < 6"^ and ap ~ 0. They are given 
as follows, 

D,\t) ^ {2M,fl'l ait,t) ^ 3 {^^y iff, 

.1.1/3 .1.2/3 

1-5i/i„((,( + p,-p)c.-(J^) OV, 

D(t,t + pa)- (^j {t' + {t + pa)^}, 

J5H^3)2^.,{t,t, -pa,pa) ~ ^ + 1^«))"^7.^(^/)"} > 

l+'m+n=l 

where all the higher order terms are non-singular. Using these expressions, Il^^i}{p) and 
(p) for < (5^ and ap ~ are expressed as 



/ 1 



where f{t,pa) is a non-singular function of t and pa. Since n^''j](p) vanishes in this limit 
and the expression of n{fj(p) has the same structure as that in continuum theory, we 
conclude that there are no (unusual) divergences such as a^p^ x oo, etc., in the vacuum 
polarization tensor. 

Finally, we investigate the logarithmic divergence. /^From the above analyses, the 
logarithmic divergence in n^,y(p) appears from the singular part in the integral for a — > 
and t ~ 0. Since nj^*2(p) is non-singular in this limit, we consider only the amplitude 
Il^^J{p). In Eq.( |4.9|) , we use the Feynman's parameter and shift the integration variables 
tfj, ^ tfj^ — ap^x. We then evaluate the contribution from the integration region < 5^, 
5 ^ 1 in the continuum limit. We then have 



2a? Jfi<5^ (27r)4 Jo {t^ + p^a^x{l - x)y r(3)2 
X ((t - paxYfiit + pa(l - x)f)hr [(/ - i)ax)-i^{if + j^a(l - x))7^] 



(4.13) 
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After some calculations, the term proportional to log a p is extracted as follows, 



Nf9 



2 



{P^^iv - P^^Pv) log a'p' . (4. 14) 



247r2 

This expression agrees with the one expected from the general analysis in Section 2. 

5 Discussion 

We have studied a perturbative aspect of a general class of Dirac operators. To avoid 
the excessive complications, we examined the simplest diagrams of the one-loop fermion 
correction to the gauge field self-energy tensor. This quantity is related to the one-loop j3 
function and also to the Weyl anomaly. We have confirmed that the perturbative analysis 
gives the correct result for any > 1 by using the relation (2.14), in accord with the 
general analysis of Weyl anomaly. This correct result is consistent with our previous 
analyses of the locality of the general operator H(2k+i) and the locality domain of \Ffj_y\ 
for H(2k+i) [12]: Also, our result does not contradict the general perturbative analysis of 
lattice theory in [26] if one remembers the locality properties of -ff(2fc+i)- We have also 
confirmed the relation (2.14) for the simplest case A; = 1 by evaluating the self-energy 
correction in terms of the operator H itself. 

When combined with the analysis of chiral anomaly [8], our present analysis shows 
that all the local anomalies are properly reproduced by our general class of operators D. 
These analyses give some confidence in the treatment of the fermionic determinant 

deti/-(deti/(2,+i))^/(2'=+^) (5.1) 

in the possible application to QCD, for example. 

At the same time, we recognized that infrared divergences may generally appear in the 
intermediate stages of perturbative calculations for finite a, which should cancel in the 
final result. This treatment of infrared divergences in perturbation theory is quite tedious 
in our generalized operator D. To avoid the infrared complications, the (nonperturba- 
tive) Wilsonian effective action, which is supposed to be free of infrared complications, is 
expected to be essential for the general operator D. As for the perturbative treatment of 
fully dynamical gauge field such as in the one-loop correction to the fermion self-energy, 
some auxiliary regulator such as the dimensional regulator may become necessary [22] for 
a reliable treatment of infrared divergences. 



A Feynman rules for the general H(^2k+i) 

We derive the Feynman rules for H{2k+\) theory (not H itself) to calculate the vacuum 
polarization at one-loop level on the basis of (2.14). //(2fe+i) has been defined by (1.6), 
(1.7) and (1.8). We expand H(2k+i) up to the second order in the coupling constant g as 
follows, 

H(2k+i) = -f^(2fe+i)o + 9H(2k+i)i + g'^H(2k+i)2 + 0{g^). (A.l) 



21 



For this purpose we first need to expand D^!''^^\x, y) in g. D^!''^^\x, y) up to the second 
order in g is given by 

DSr'^'\x, y) = ^" [Xo{p)5p{p -q)+ X,{p, q) + X,{p, q) 

, (A.2) 



where 



2k+l \2fe+l /^x2fe+l 



-7^sinp^aj + ( - ^(1 - cosp^a) 1 - 



a J 



X,ip,q) ^ r 7PrMP-^-<i)9A,{t)X,,ip,q,t), (A.4) 
xX2^,u{p,q,ki,k2). 

(A.5) 

Here we defined 

Xi^i{p, q,t)^ (-^p sinppo) (i'j^ cos(g + ^)^a) (-7^ sin qpo) 

l+m=2k I ^ ^ Z / \a / 

+ (^^51(1 -COS ppa)j (^rsin(5 + ^)^aj |^^^(l-cosgpa)j l, (A.6) 



and 



Z+m+n=2fe-l I ^ \ ^ / 

X f -7p sin(p - ki)pa] ( ^7^ cos(? + —)^a j f -7p sin g^a j 
+ (^^ - cosppa) j sin(p - y 

X ^(1 - cos(p - ki)pa)j (^r sin(g + y (^^ $^(1 - cos 5^0)^ | 
+\ H |iQ7pSinppa^ f-ia7^5^^sin(g + ■^i^-^)^a'j Q7pSingpa^ 
+ ^^^(1 - cosppa)j ^ar^^,. cos(g + ^^-^-^)^aj ^(1 - cos gpa)j 



(A.7) 
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where 



2fc+l 



M(p) 



^(1 - cosppo) 



/mo 
V a 



2k+l 



(A.8) 



and /, m, n are the nonnegative integers and 5p is the periodic lattice delta function. A^{k) 
is defined by 



A^{k)e 



L.(2vr)4 " 

and has the properties Aj^{k) = A^{-k), A^{k + = {-lY^'A^{k) {I : integer). 



(A.9) 



Next we want to expand the factor 1/ \J {D^^'^^y {d'"^'^^^) in -ff(2fc+i) ( |1.6|) . But it 



is very complicated to perform the weak coupling expansion of 1/ y {D^^^^)'^{D^^^^) 
directly, if not impossible. Therefore we use the following identity: 



Vxtx 



dt 



oo TT t2 + XtX ' 



(A. 10) 



The weak coupling expansion of the integrand on the right-hand side can be readily 
performed. After some calculations, we obtain 



(2k+l) 



d'^p f d'^q 



(27r)4 7 (27r)4 

uj{p)+uj{q)} I 

uj{p)+uj{q) J \ 
1 1 



^ipx-iqy 



,1 

2 



1 + 



uj{p) _ 



6p{p - q) 



Xi [p, q) -— (p, g) — — 

uj{p) uj{q) 

^2 (p, g) VP, q) 



uj{p) 



uj{q) 



t [uj{p) + iu{q) } [uj{p) + iu{t) j [Lu{t) + Lu{q) j 
x{-Xo{p)Xl{p,t)X^{t,q) 

-Xi(p, t)Xl{t)X,{t, q) - X,{p, t)Xl{t, q)Xoiq) 

uj{p)uj{t)uj{q) 



(A.ll) 



where 



uj{p) 



\ 



-sm 



2fc+l 



+ 



2fc+l 



COSPpflj 



mo 
a 



2fc+l 



(A.12) 
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and = 75X75. We write H(^2k+i)o{x,y), H^2k+i)i{x,y) and H^2k+i)2{x,y) as follows, 

H^2k+i)iix, y) = a'Y.[ e'^^-''^y5p{p - q - t)A,{t)H^2k+i)i,{p. q, t), (A.14) 

H(2k+i)2ix, y) = a^Y. I e'P''-''^y6p{p - q - h - k2)A^{h)A^{k2)H(^2k+i)2,,u{p, q, h, ^^2), 

(A.15) 

where = Then from Eq.( |A.ll|) , H^2k+i)o{p), H^2k+i)if,{p,q,t) and 

H{2k+i)2fiu{p, q, h, h) can be written as 



1 



H(2k+i)i,i{p,q,t) 



t75 



2'^a.(p)+u;(g) 



xjx,,(p,g,t)-^Xt,(p,g,t)-^°('^) 



uj{q) 



1 1 

H(2k+i)2f,u{p, q, ki, ^2) = 7t75 . X , ^ 

2 c<j(p) + u!{q) 



X 



X2^u{p,q,ki,k2 



^^4^X1 {p,q,ki,k2] 



Xojq) 
uj{q) 



(A.16) 



(A.17) 



+ iu{p - ki)){uj{p - ki) + 
X {Ari^(P,P - ki, ki)Xl^{q + k2,q, k2)Xo{q) 
+Xif,{p,p- ki, ki)Xl{p - ki)Xiy{q + k2,q, k2) 
+Xo{p)Xl^{p,p - ki, ki)Xiy{q + k2,q, ^2) 
iu{p) + uj{p - ki) + ijj{q) 
uj{p)uj{p - ki)uj{q) 
xXo{p)Xl^{p,p- ki, ki)Xo{p - ki)Xl^{q + ^2, g, ^^2)^0 

(A.18) 

Note that q = p - t in H(2k+i)i,,{p, q, t) and q = p - ki - fcs in H(2k+i)2,,v{p, q, ki, ^2). 

/^From this weak coupling expansion we can derive the Feynman rules for H{2k+i) 
theory, which are necessary for the one-loop analyses. To make the structure of the 
divergences in one-loop amplitudes explicit we derive the Feynman rules for ^2k+i H(2k+i) = 
(75^)^'^+^. Using Eq. (|A.16|) , the fermion propagator D^^{p), where Dq^{p) stands for the 
inverse of free {'j^D)'^^'^^ , is written as 



^w{p) + M{p) 



75 



(A.19) 
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where = Z]/^ 7/^ sin ap^ and = Z^^sin^ap^. Using Eq. (|A.17| ), we assign the following 
factor to the fermion-gauge field three-point vertex depicted in Fig. 2 



9 



2a2fc+l 



1 



uj{p)+uj{q) 



Xi,(p,g,t)-^Xt,(p,g,t)-^ 



(A.20) 



where are SU{N) generators. Using Eq.( |A.18D , we assign the following factor to the 
fermion-gauge field four-point vertex depicted in Fig. 3 



9' 



uj{p)+uj{q) 



X2i,u(,P, q, ki, - (p, q, ki, /cs)^"*^^'' 



uj{p) 



uj{q) 



{uj{p) + uj{p - ki)){uj{p - ki) + uj{q)) 
X {Xi^(p,p - ki, ki)Xl^{q + k2,q, k2)Xo{q) 
+Xi^{p,p- ki, ki)Xl{p - ki)Xi^{q + k2,q, /C2) 
+Xq{p)XI^{p,P - ki, ki)Xi^{q + k2,q, k2) 

uj{p) + uj{p - ki) + uj{q) 



uj{p)uj{p - ki)uj{q) 
xXo{p)Xl^{p,p- ki, ki)Xo{p - ki)Xl^{q + fca, q, k2)Xo{q)}] 
+ {A,fi,ki^ B,u,k2) (A.21) 
where we have imposed the Bose symmetry for gauge fields. 



p,a 



q,b 



Figure 2: Fermion-gauge field three-point vertex 



In Section 3 of the present paper, we calculate the vacuum polarization tensor at 
one-loop level by using these Feynman rules. 



B Feynman rules for the operator H with k=l 

In this Appendix we derive the Feynman rules for the operator H to calculate the vacuum 
polarization at one-loop level. For simplicity we consider the case with k = 1. For the 
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A,iJ,,ki B,p,k2 



P>a q,b 
Figure 3: Fermion- gauge field four-point vertex 

Feynman rules for the operator i?(3), we refer to Appendix A. We expand H and i/(3) up 
to the second order in the coupling constant g as follows, 

H = Ho + gH,+g''H2 + 0{g^) (B.l) 

D = --f5Ho + --f59Hi + -^59^H2 + 0{g^) 
a a a 

Hi3) = i/(3)0 + + /^i'(3)2 + 

= {Ho + gH^ + g^H2 + Oig^)y H^) 
= Hl + g Y: KH.Hl-^ 

0<m<2 

\0<m<2 0<Z+m<l 

(B.2) 

Now we want to know the expressions of Hq, Hi, H2. We have obtained Hq in our previous 
paper[12][9]. In momentum space Hq is written as 

2 2 

Hoip) = 75Q)'(^)'|(a;(p) + M(p))i-(a;(p)-M(p))i^|, 

= -(24))'(otr)'{-@^^-^^^^"W' 

(B.3) 

where ^ = 7/i sin ap^u and Sp = S^^sin^ap^. One can easily check that i7(3)o(p) = 
Next we derive the expression of Hi. We write Hi as follows, 

//i = a^E/ e^"--<'^<Jp(p-g-t)A^(t)i/i^(p,g,t), (B.4) 
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where J^ g ^ = J ^ J ^ J Using Eq.(|B3) at first order in g, we obtain 

0<m<2 

= D{p, q)Hi^{p, q, t) + Ho{p)H^^{p, q, t)Ho{q), 

(B.5) 

where 

D{p,q) = Hoipf + Hoiqf, (B.6) 



'1 A I Mip)\ 



Now we consider an ansatz as 

Hl^,{p, q,t) = - (d{p, g)i/(3)i;.(p, g, t) - Ho{p)H^3)i^{p, q, t)i/o(g)) • 

(B.8) 

Substituting this ansatz for Hi^ into Eq.( |B.5| ), we easily obtain the expression of a as 

aip,q) = D\p,q)-Hl{p)Hl{q) 

= J2 H,{pf"^H,{qf^^~"^\ (B.9) 

0<TO<2 

Thus we have obtained the expression of Hi. 

Further we derive the expression of H2, performing the similar procedure as deriving 
Hi. We first write H2 as follows, 

H2 = e'^--''^y5p{p-q-ki-h)A^iki)A,{h)H2,.{p,q,ki,k2). 

(B.IO) 

Using Eq.( p.2|) at the second order in ^f, we obtain 

Hc,)2AP^q,ki,k2)= E H^{prH2,AP.q.kiM)Ho{qf-"' 

0<m<2 

+ Yl Ho{p)'Hi,{p,p- ki, ki)Ho{p - kirnUp - ki, q, k2)Ho{qy~'-"' 

0</+m<l 

= D{p,q)H2^„y{p,q, /ci, ^2) 

+Ho{p)H2^u{p, q, ki, k2)Ho{q) 

+Hif,{p,p - ki,ki)Hi^{p - ki,q, k2)Ho{q) 

+Hi^{p,p- ki,ki)Ho{p- ki)Hi^{p- ki,q, ^2) 

+Ho{p)Hi^{p,p- ki, ki)Hi^{p - ki, q, /C2) (B.ll) 
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Now we define 

— {Hifj_Hii,Ho + Hi^HqHii, + H^Hi^Hiy} . 

(B.12) 

Using -ff(3)2^,y) Eq.( |B.ll| ) is written as 

i/('3)2^, = DH2^. + Hoip)H2,Mq)- (B.13) 
The structure of this equation is the same as that of Eq.( |B.5| ). Therefore H2fj,u is written 



H2f^u = ^ J {Dip, q)H[3)2i,u - Ho{p)H[3)2f,uHo{q)) ■ 



(B.14) 



Finally we obtain the expression of -^2/^1/ in momentum space as follows, 

H2^u{p,q,h,k2) = ■^^^^^—:^[D{p,q)H(^3)2f,u{p,q,ki,k2) 

-Ho{p)H(^s)2f,u{p, q, ki, k2)Ho{q) 
-Ho{qyHi^{p,p - ki, ki)Hi^{p - ki, q, k2)Ho{q) 
-D{p,q)Hif,{p,p - ki, ki)Ho{p - ki)Hi^{p - ki,q, /C2) 
-Ho{pyHo{p)Hi^{p,p - ki,ki)Hi^{p - ki,q, ^2) 
+Hq{j))Hi^{p,P - ki, ki)Ho{p - ki)Hi^{p - ki, q, ^2)^^0(5)} • 

(B.15) 

/^From this weak coupling expansion we can derive the Feynman rules for D in the case 
oi k = 1, which are necessary for the one-loop analysis. Using Eq.( |B.3| ), the fermion 
propagator Dq^{p) is written as 

Using Eq.( p.8| ) and Eq. ( P.9|) , we assign the following expression to the fermion-gauge field 
three-point vertex depicted in Fig.2, 

9rpA 1 



X 



a '"'a{p,q) 

\D{p,q)lbH(^i)i^,{p,q,t) --fr^Ho{p) (75^(3)iM(P,9,^))S5^o(g)| , 

(B.17) 
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where are SU{N) generators. Using Eq.( [B.15|) , we assign the following expression to 
the fermion-gauge four-point vertex depicted in Fig.3, 

{T^T^)ba— r {d{p, g)75if(3)2^^(p, g, /ci, k2) 

a a[p, q) ^ 

X'y5Hif,{p,p- ki,ki) {j^Hi^{p-ki,q, k2))^ -f^Hoiq) 
-D{p,q) 

x-f5Hif,{p,p - ki,ki) {j^Hoip - ki)y 'j^Hi^ip - ki,q, /ca) 
-Hoipf 

x^bHoip) {-f^Hi^i^p.p - ki.ki))^ -i^Hi^{p - ki,q, /cg) 
+l5Hoip) {-f5Hif,{p,p- ki, ki))^ -i^Hq{p - ki) 

X (75^i.(p - ku q, fcs))^ 75^o(g)} + {A /i, ki ^ B, u, k2). (B.18) 
We perform a one-loop calculation in Section 4 on the basis of these Feynman rules. 
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